Coherence properties of a 2D trapped Bose gas around the superfluid transition 
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We measure the momentum distribution of a 2D trapped Bose gas and observe the increase of the 
range of coherence around the Berezinskii-Kosterlitz-Thouless (BKT) transition. We quantitatively 
compare our observed profiles to both a Hartee-Fock mean-field theory and to quantum Monte-Carlo 
simulations. In the normal phase, we already observe a sharpening of the momentum distribution. 
This behavior is partially captured in a mean-field approach, in contrast to the physics of the BKT 
transition. 
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The characteristics of a Bose gas in the degenerate 
regime are greatly dependant on the dimensionality of 
the system. In contrast to its three-dimensional (3D) 
counterpart, a two-dimensional (2D) interacting Bose gas 
does not present true long range order at low tempera- 
ture [l|, Q . Instead, it undergoes a superfluid-to-normal 
fluid transition described by the BKT theory 0, 3] ■ This 
superfluid transition has been observed in *He films 
Josephson coupled superconducting arrays and re- 
cently in dilute ultracold atomic gases 0. 

In ultracold 2D trapped Bose gases, the apparition of 
superfluidity was predicted to be preceded by the reduc- 
tion of density fluctuations 0, Q and by the increase of 
coherence , [lH . Experimentally, great efforts were de- 
voted to the study of the density distribution n(r) of 2D 
gases, which gives direct access to the phase space den- 
sity and to the density fluctuations. It has permitted the 
study of the scale invariance [l2| and universality in 2D 
Bose gases as well as its thermodynamics [iJI- The 
study of the first order correlation function gi is possible 
in real space through interferences between two clouds of 
atoms [71, [15 1 which revealed its expected algebraic de- 



cay in the superfluid regime 0. The Fourier transform 
of gi is also naturally embedded in the momentum dis- 
tribution 1^ ill- For example, the narrow peak in the 
density distribution after time of flight is frequently used 
as indicator of the slow algebraic decay of the coherence 
in the superfluid regime [H, [H, [l^ • 

In this paper, we present a detailed analysis of the 
experimental momentum distribution of 2D Bose gases 
closer to a genuine 2D gas than previous studies [I3 • In 
addition, we fit our profiles using a Hartree-Fock mean- 
field (HFMF) model already used in density space 
but extended here for the momentum distribution, and 
compare our results to quantum Monte-Carlo (QMC) 
simulations . As we increase the phase space den- 

sity, we observe a progressive narrowing of the momen- 
tum distribution and a change of shape from an approx- 
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FIG. 1. Scheme of the experimental setup. The atom cloud, 
initially trapped in a harmonic 2D trap with frequencies 8 Hz 
X 15 Hz X 1.5 kHz, is dropped for 83.5 ms before shining 
a fluorescence beam situated 3.4 cm below the initial posi- 
tion. The fluorescence signal is recorded from above with an 
electron-multiplying-CCD camera. 



imately gaussian to a peaked distribution. In particular, 
at the superfluid phase transition, the momentum distri- 
bution is already peaked and its width reduced by a fac- 
tor ^ 5 compared to a Boltzmann thermal gas. Whereas 
this behavior is partially captured in a mean-field ap- 
proach where the 2D Bose statistics plays a central role, 
correlation effects beyond mean field are cleary visible in 
the full momentum profile. 

Experimentally, the 2D Bose gases are prepared as fol- 
lows. We prepare a 3D cloud of ^^Rb atoms cooled by 
all-optical runaway evaporation in a crossed dipole trap, 
as described in |2l|. During the evaporation process, we 
turn on a uniform magnetic gradient to select a single 
hyperfine state (|F = 1, mp = 0)) [l^l- We then trans- 
fer the cloud to a 2D trap. For that, we adiabatically 
ramp on a blue detuned laser at 767 nm in a TEMoi-like 
mode, providing an intensity node in the horizontal plane 
of the atoms (see FigllJ. It induces a strong confinement 
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in the vertical direction, with an oscihation frequency of 
LOzl^T^ = 1-5 kHz, corresponding to a dimensionless inter- 
action strength g = gm/h? = -v/STras/a^ = 0.096, where 
fls = 5.3 nm is the 3D scattering length, az = ^h/miUz 
is the harmonic oscillator length in the vertical direction, 
h is the reduced Planck constant, and m is the atomic 
mass. 

Horizontally, the 2D trap is made of a laser beam at 
1565 nm with a waist w = 200 /zm and tilted by ~ 30° 
with respect to the horizontal plane (see FiglT]). Its os- 
cillation frequencies are uj^/'^t^ = 8 Hz, ujy/2'K = 15 Hz. 
The atom number N is varied from 2 x 10^ to 6 x 10^ by 
changing the number of atoms initially loaded. The fi- 
nal temperature remains approximately constant at T = 
64.5 ± 2.0 nK. The two dimensional character of our 
experiment is given by k^T /fnoz = 0.90 where fee is the 
Boltzmann constant, which results in having ~ 70 % of 
the atoms in the ground state of the vertical harmonic 
oscillator. 

We let the atoms thermalize for 500 ms in the final 
trap before probing the momentum distribution function 
through time of flight imaging (TOF). In the first mil- 
liseconds of expansion, the gas expands predominantly 
in the vertical direction and the interaction energy is 
thus quickly released in this direction. As a result, the 
horizontal momentum distribution (in the xy plane) re- 
mains unchanged during the expansion. After an expan- 
sion time of ixOF = 83.5 ms much longer than l/ux and 
\/u)y^ the horizontal density distribution reflects the ini- 
tial momentum distribution [2^. Finally, two circularly 
polarized saturating retro-reflected laser beams, resonant 
with the |5Si/2, F = 2) -> 
and with the |5Si/2, F = l) 
tion allow the atoms to fluoresce for 100 ^s. The fluores- 
cence signal is recorded on an electron-multiplying-CCD 
camera placed along the vertical axis, thus imaging the 
horizontal profile. 

To analyze the experimental data, we perform an az- 
imuthal averaging of the single pictures, as described in 
[l3 |. in order to extract radial profiles corresponding to 
the momentum distributions n(|k|), where k is the atom 
wavevector. Experimentally, we weight every point of 
the profile by its experimental standard deviation (com- 
ing from the photon shot noise, the camera dark noise 
and the number of averaging points). Such profiles are 
presented in Fig.[2] For a low atom number (see Fig. [2^), 
the distribution is approximately gaussian and relatively 
broad. For a higher atom number (Fig.[2j;-d), the curve 
progressively peaks and deviates from a gaussian. At 
high atom number (Fig.[2^-b), a sharp feature devel- 
ops at very low momentum and grows with increasing 
atom number. Its width tends to a constant and agrees 
with our resolution in momentum space (half width of 
~ 0.3/im^^) coming essentially from our imaging reso- 
lution and from the initial size of the cloud. Since the 
peak at low momentum develops on a distribution that is 
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FIG. 2. Radial profile of the momentum distribution for 
five different atom numbers at constant temperature T = 
64.5 ± 2.0 nK. The atom number A'^ is given in unit of the 
critical number for the ideal gas Bose-Einstein condensation 
A^c ~ 3 X 10* and the superfluid phase transition is expected 
for N/Nc ~ 1.18 (see text). In each plot, we present the exper- 
imental data (black dots), the mean- field profiles (blue dashed 
line) resulting from the fit in the wing (i.e. for k > 2/im~^) 
and the corresponding quantum Monte-Carlo profiles (green 
continuous line) with the same temperature and the measured 
atom number. 



not a gaussian, it is hard to precisely point when it first 
appears. 

In order to attribute a temperature to each profile, 
we have to rely on a model. More precisely, we fit the 
Hartree-Fock mean-field model to the wings of the mo- 
mentum distribution data. This model has already been 
used for in-situ density profiles [13, 12 , 2^ but needs to be 



extended to get the momentum distribution. We proceed 
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as follows. The density distribution in HFMF theory, in 
the local density approximation, reads 
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(1) 



where /3 = l/ksT, fi{r) = /xq— "T-w^a;^/2— is the 
local chemical potential with fiQ the chemical potential 
at the trap center. After integration, this leads to the 
following equation relating fi{r) and n{r) [2^ : 



(2) 



where AdB = ^j2'Kh? /mk^T is the thermal de Broglie 
wavelength. Writing the momentum distribution: 
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and making the change in variables from (x, y) to ^, we 
directly obtain n(k) after integration. By taking into 
account the thermally populated vertical levels but ne- 
glecting the interaction in these levels, which is justified 
since the densities in the excited levels are small, we fi- 
nally calculate the function : 
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(4) 

Where is determined numerically by solving Eq.([2]), 
and where oj^ = a;^ + a;^. For given values of /^o and T, 
we can calculate JT-dkl) and then use this function to fit 
to the experimental data. 

More precisely, we fit the wings of the data where 
we expect beyond mean-field effects to play little role 
and we can extract temperature and chemical potential. 
This requires an accurate calibration of the atom detec- 
tion efficiency, which is a rather difficult task. We have 
performed QMC simulations based on numerically exact 
path-integral algorithm [2^ [2^ to calculate the mo- 
mentum profiles, so that we can adjust the calibration of 
the atom detector. Within the experimental error bars, 
it agrees with an independent calibration using the 3D 
condensation threshold. Moreover, by fitting a mean- 
field model to the wings of QMC simulations, we find 
the accurate temperature and thus validate our fitting 
method. 

In FigI21 we plot with each profile the fitted mean- 
field distribution (blue dashed line). We observe that 
the mean-field fit reproduces well the experimental data 
at low atom number, while it fails accounting for the 
central part of the profiles when the number of atoms 
increases. However, the QMC simulations for the fitted 
temperatures and the experimental atom numbers (green 
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FIG. 3. a. Half width at half maximum of the profiles, nor- 
malized to the HWHM of a gaussian profile for the same 
temperature, plotted versus the number of atoms normal- 
ized to A^c- b. Fraction of atoms in the central pixel of the 
image No/N. Blue dashed line : mean-field prediction for 
T = 64.5 nK. Green continuous line: Monte-Carlo data for 
T = 64.5 ± 0.3 nK. Red vertical dot-dashed line : superfluid 
transition. 



continuous lines in FiglJ) are in agreement with the ex- 
periment in all regimes. For both the calculated mean- 
field and QMC profiles, we take into account the finite 
resolution of our imaging system. 

In order to quantify the degeneracy of the gas, we cal- 
culate N/Nc, where 



(5) 



is the critical atom number for a non-interacting Bose gas 
for our trap parameters, and Q2 (x) = /n^. From 

classical field calculations 1131 j the BKT transition 
is expected at a central density of the lowest vertical level 
ln(380.3/o)/Ajg with small quantum corrections of order 
g/'^dB [lOl, ISOf. Integratin g th e dominating classical field 
corrections to mean-field |10| yields (A/A^c)bkt ~ 1-18 
as critical value of the particle number. Close to this 
value, the experimental profiles are already peaked at 
low momentum (Figl^b-c). 

In order to analyze our experimental findings further, 
we now consider two quantities : the half width at half 
maximum (HWHM) and the fraction of atoms in the cen- 
tral pixel Nq/N. The advantage of those two quantities is 
that they are model-independent and quantify the degree 
of coherence of the gas. The HWHM gives a measure of 
the inverse of the coherence length whereas the fraction 
of atoms in the central pixel is related to the fraction of 
atoms which are coherent on a length scale larger than 
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^ 5 /iiTL. In FiglSl wc plot these two quantities as a 
function of N/Nc- The HWHM is normahzcd to the one 
of a gaussian distribution : 2y^n ln(2)/AdB = 4.0 ^m^^ . 
Aheady for N/N^ « 0.5, the normahzed width of the 
momentum distribution starts to decrease from 1. In 
other word, the distribution is not a gaussian any more 
and tends to peak at low momentum. At the superfluid 
transition point {N/Nc « 1.18), the HWHM has already 
decreased by a factor of 5. After the transition, the 
HWHM saturates to a value corresponding to the imag- 
ing resolution. 

The change of shape in the momentum distribution 
is also reflected in the fraction of atoms in the central 
pixel (Fig. 3b). This value increases smoothly but quite 
dramatically as the number of atoms increases. For low 
atom number, it corresponds to the decrease of the width 
and for high atom number, while the HWHM saturates, 
the fraction of atoms in the central pixel keeps growing, 
reflecting the increasing coherence of the gas with atom 
number. It is interesting to note that the superfluid tran- 
sition does not appear as a sharp feature in our data on 
the coherence properties of the gas, unlike what is ex- 
pected for the superfluid fraction 
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In addition to the experimental points, we also plot 
Monte-Carlo simulations and mean-field calculations for 
our experimental conditions. The Monte-Carlo simula- 
tions show a good agreement with our experimental find- 
ings. The mean-field results coincides with the Monte- 
Carlo simulations at low atom number up to about 
N/Nc ~ 1- For high degeneracy parameters N/Nc > 1, 
the mean-field model underestimates the height of the 
coherence peak showing that beyond mean-field effects 
become important. It is remarkable that the mean-field 
approximation captures the initial increase of the coher- 
ence length. This effect is thus not in direct link with the 
physics of the BKT phase transition for which the theory 
is inherently beyond mean- field. In fact, even in a non- 
interacting trapped Bose gas, for which calculations are 
exact, there is also an increase of the coherence length 
before the Bose-Einstein phase transition. 

The presence of a low momentum coherence peak can- 
not be considered as a signature of the BKT phase tran- 
sition. An accurate signature however is the slow al- 
gebraic decay of the first order coherence function (as 
l/r" with a < 0.25) Q- This decay can in principle be 
observed in the momentum distribution and indeed our 
QMC calculations show a change of slope (to negative) 
in the function k'^~°'n{k) at low momentum, k < Ads for 
N > 1.27Vc, together with the onset of superfiuidity. Un- 
fortunately, our experimental resolution is not sufficient 
for direct observation of this feature. 

In conclusion, we have studied the momentum distri- 
bution of a trapped interacting 2D Bose gas. Our finding 
are in qualitative agreement with previously published re- 
sults but we analyze our profiles in more de- 
tails and compare them to mean-field and Monte-Carlo 



calculations. In particular, we show that the momen- 
tum distribution narrows progressively and well before 
the BKT phase transition. Our detailed characteriza- 
tion of the 2D trapped Bose gas momentum distribution 
will be a useful tool for further studies of the remarkable 
properties of 2D gases of ultracold atoms. 
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